Initialization and Readout of Spin Chains for Quantum Information Transport 
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Linear chains of spins acting as quantum wires are a promising approach to achieve scalable 
quantum information processors. Nuclear spins in apatite crystals provide an ideal test-bed for the 
experimental study of quantum information transport, as they closely emulate a one-dimensional 
spin chain. Nuclear Magnetic Resonance techniques can be used to drive the spin chain dynamics 
and probe the accompanying transport mechanisms. Here we demonstrate initialization and read- 
out capabilities in these spin chains, even in the absence of single-spin addressability. These control 
schemes enable preparing desired states for quantum information transport and probing their evo- 
lution under the transport Hamiltonian. We further optimize the control schemes by a detailed 
analysis of 19 F NMR lineshape. 

PACS numbers: 03.67.Hk, 03.67.Lx, 75.10.Pq, 76.90.+d 



I. INTRODUCTION 

Control over small quantum systems and the ability to 
perform simple quantum algorithms have been demon- 
strated on a variety of physical systems ranging from 
trapped ions [1] and electrons [2], to neutral atoms and 
molecules in optical lattices [3], to superconducting cir- 
cuits [4] and semiconductor quantum dots [5] , to nuclear 
and electronic spins [6-8]. Although algorithms involv- 
ing more than one qubit have been executed [9-14], a 
vital requirement for a quantum computer -scalability 
while preserving fidelity- has not yet been achieved in 
any physical system. The use of linear chains of spins as 
quantum wires to couple basic memory units is a promis- 
ing approach to address this issue [15, 16]. These spin 
chains have the ability to transmit quantum information 
via the free evolution of the spins under their mutual 
interaction [17-22]. While advances in fabrication tech- 
niques have made physical implementation of spin wires 
possible [23-26], the level of precision available is not 
yet adequate. Therefore, natural systems such as crys- 
tals where spin position is precisely set by nature are a 
preferred choice for exploring such applications. 

Owing to their unique geometry, nuclear spin systems 
in apatite crystals have emerged as a rich test-bed to 
probe quasi-one-dimcnsional (ID) dynamics, including 
transport and decoherence [27-30]. The crystal wherein 
19 F (or 1 H) nuclei are aligned along one axis, emulates 
a collection of ID chains. The dynamics of these spin 
chains have been studied by various nuclear magnetic 
resonance (NMR) techniques [31-33]. In our previous 
work [20, 34], we have shown that the natural dipolar 
interaction among the spins can be manipulated via the 
available collective control to simulate the Hamiltonian 
driving quantum transport. The lack of single-spin ad- 
dressability in this ensemble systems however would seem 
to prevent creating and measuring a single-spin excita- 
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tion as required to study transport. Still, we demon- 
strated experimentally [27, 28] in the Fluorapatite (FAp) 
system that one can prepare the spin system in an initial 
state in which polarization is localized at the ends of the 
spin chain, a state that well simulates the conditions for 
spin-excitation transport [34]. 

In this paper we take further steps toward enabling 
the experimental study of quantum transport in spin 
chains: we introduce an experimental technique to read 
out the spins at the chain extremities and we show 
how to prepare a two-spin encoded state that is able 
to transfer quantum information. We use these initial- 
ization and readout techniques to study the dynamics 
under the transport Hamiltonian (the Double Quantum 
(DQ) Hamiltonian [30, 35]). Additionally, we probe 
the spin chain dynamics by creation and evolution of 
multiple quantum coherences [36], which present well- 
characterized state dependent signatures. We use both 
these techniques to demonstrate preparation as well as 
readout of spins at the chain ends, thus verifying impor- 
tant required tasks towards simulation of quantum infor- 
mation transfer. We further validate the addressability of 
ends of the chains by a detailed analysis of 19 F lineshape 
in FAp. 

These techniques will make it possible to explore er- 
rors affecting the transport fidelities as well as control 
schemes to mitigate them, in an experimental setting, 
where the interactions among spins are not limited to 
the ones tractable by solvable models. 



II. TRANSPORT IN MIXED-STATE SPIN 
CHAINS 

A. Spin chain dynamics 

Linear chains of spin- 1/2 particles have been proposed 
as quantum wires to transport quantum information be- 
tween distant nodes of a distributed quantum computing 
architecture. Transport can occur even in the absence 
of individual control of the chain spins, as it is medi- 
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ated by the spins mutual interactions. In the most stud- 
ied model, energy conserving spin flip-flops (mediated by 
the isotropic XY Hamiltonian) drive the transport of a 
single-spin excitation [17-22]. This model has been re- 
cently extended to the case where the initial state of the 
chain cannot be perfectly controlled, and thus it is found 
in a mixed-state rather than its ground state [34, 37-40] . 

Spin chains that are in a maximally mixed-state are 
particularly interesting from the point of view of exper- 
imental study of quantum information transfer. This 
state, corresponding to the infinite temperature, can be 
easily achieved experimentally and has been shown to 
provide a direct simulation of pure state transport [34]. 
Additionally, extension to mixed-state chain enables us- 
ing the so called double quantum (DQ) Hamiltonian: 

to drive transport, although it does not conserve the 
number of spin excitations. This Hamiltonian can be 
easily obtained from the natural dipolar Hamiltonian 
with only collective control [20, 34] and it is related to 
the isotropic XY Hamiltonian (which instead cannot be 
generated from the dipolar interaction) via a similarity 
transformation. The extension to mixed-states and to 
DQ Hamiltonian open the possibility to study experi- 
mentally quantum information transport in nuclear spin 
chains with NMR techniques. Under our experimental 
conditions (strong external magnetic field, Bo=7T and 
room temperature), the initial equilibrium state is the 
Zeeman thermal state, 

,4(0) cx exp(-eo*)«]L-eff*, (2) 

where a z = J2j a j an d e = iBo/ksT. Since the identity 
does not evolve and does not contribute to the signal, we 
will focus on the deviation of the density operator from 
the maximally mixed state, 5p ~ p — 1. In the absence 
of individual spin addressing, transport within a chain 
can be studied by preparing a polarization excess at one 
end of the chain, that is, a state where one spin at the 
chain extremities is polarized while the remaining spins 
are fully mixed, Sp ~ u\ . Because of the symmetry be- 
tween the two chain ends, the state we can prepare exper- 
imentally [27, 28] , which we call "end-polarized state" , is 
given by 

<Wi(0) = af+a z N . (3) 

The end-polarized state simulates the dynamics of a 
single-spin excitation in a pure-state spin chain. This 
state can transfer a bit of classical information by en- 
coding it in the sign of the polarization. This encoding 
is, however, not enough to transfer quantum informa- 
tion, which requires additional transfer of information 
about the phase coherence of a state. A two-qubit en- 
coding allows for the transport of a bit of quantum in- 
formation [34, 39]. For transport via the DQ Hamilto- 
nian, this encoding is given by the basis |0)^ q = |00) and 



|l) L q = 1 1 1) . The operator basis for transport via mixed 
states under DQ Hamiltonian is thus given by: 

dq _ dq _ 



Starting from any of the above initial states, the evo- 
lution under the DQ Hamiltonian directly simulates the 
transport dynamics within a chain. 

In the limit of nearest neighbor (NN) coupling only, 
the evolution under DQ Hamiltonian is exactly solvable 
by invoking a Jordan- Wigner mapping onto a system of 
free fermions [30, 41, 42]. The resulting dynamics of the 
observables we analyzed in our experiments have been 
reported in the literature (see e.g. [30, 34]) and are re- 
viewed in Appendix A for completeness. Isolated, lin- 
ear spin chains with NN couplings is an accurate model 
for the experiments, given the experimental timcscales 
used [28]. Comparison of the theoretical model with the 
experimental results thus allows us to validate our ini- 
tialization and readout methods. 

To gather more insights into the states generated by 
the evolution under the DQ Hamiltonian, we experimen- 
tally measured multi-spin correlations via multiple quan- 
tum NMR experiments. This experimental technique 
(see Appendix B) allows measurement of multi-spin state 
dynamics by indirectly encoding their signatures into sin- 
gle spin terms. The dynamics of quantum coherence in- 
tensities can as well be calculated analytically in the limit 
of NN couplings [20, 34, 42] and we review these results 
in Appendix B. 

B. Preparing and reading out desired spin states 

To probe the quantum transport dynamics it is nec- 
essary to prepare the spins at the ends of the chain in 
a perturbed state (such as Sp cn d), which is then left to 
evolve under the transport Hamiltonian. Furthermore, 
to observe the transport, measurement of the end chain 
spins would be desirable. Unfortunately, in a system of 
dipolarly coupled homonuclcar spins (such as FAp) it is 
not possible to achieve frequency or spatial addressability 
of individual spins. Still, here we show that we can ap- 
proximate these preparation and readout steps with the 
control at hand, combining unitary and incoherent spin 
manipulation. The key observation is that even in the 
absence of frequency addressability, the dynamics of the 
end-chain spins under the internal dipolar Hamiltonian is 
different from the rest of the spins in the chain. Indeed, 
the spins at the ends of the chain are coupled to only one 
nearest neighbor whereas spins in the rest of the chain 
have two nearest neighbors. This fact can be exploited 
to experimentally prepare the spins at the chain ends in 
a desired state [27, 28] as well as subsequently read out 
this state as explained below. 

When the initial thermal equilibrium state is rotated 
to the transverse plane by a it/2 pulse, we create a state 
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X^fcLi <t* which evolves under the internal dipolar Hamil- 
tonian. Due to fewer numbers of couplings with neigh- 
boring spins, the spins at the end of the chain have slower 
dynamics as compared to the rest of the chain. Thus, one 
can select a particular time t\ such that whereas the state 
of the spins at the ends is still mainly a x , the rest of the 
spins have evolved to many-body correlations. A second 
7r/2 pulse brings the magnetization of the end spins back 
to the longitudinal axis while an appropriate phase cy- 
cling scheme cancels out other terms in the state. We 
used the following pulse sequence and appropriate phase 
cycling scheme to select the ends of the chains, 



tt/2| q 



ti 



tt/2| 



(PI) 



with a =x,y to average out terms that do not com- 
mute with the total z-magnetization. For FAp crystals 
we found that the optimal t\ time (which we will call 
"end-selection time") is given by 30.3 ps [27, 28] . Fur- 
ther details on how we optimized this time are given in 
Section HID. 

The end-selection scheme forms the basis for a strat- 
egy to prepare other states, presented in Eq. 4, required 
for quantum information transport. In order to prepare 
these encoded states experimentally, we use the follow- 
ing scheme. We first prepare the end polarized state 
<5pend(0) = <j\ + g z n and then let the system evolve un- 
der DQ Hamiltonian for a short time. Applying a double 
quantum filter by an appropriate phase cycling scheme 



T !)L 



J !,L 



selects the desired state: Sp y L 

Similarly, a 7r/4 collective rotation around z, prior to the 
double-quantum filter, is needed to select the Sp^. opera- 
tor. 

A similar combination of unitary and incoherent spin 
control can be used to read out the spins at the end of 
the chain. In inductively measured NMR, the observable 
is the collective magnetization of the spin ensemble. To 
simulate the readout of a different observable, the de- 
sired state must be prepared prior to acquisition. The 
sequence PI can be used for this purpose: we call this 
the "end readout step" . We note that the sequence used 
for readout cannot in general be a simple inversion of the 
end-selection step (which is not a unitary -reversible- 
operation). A proper phase cycling should ensure that 
the state prior to the end-selection sequence has con- 
tributions mainly from population terms (oc er^). This 
property is verified for the states produced from evolu- 
tion under the DQ Hamiltonian in ID systems, making 
the PI sequence effective for the end- readout step; a more 
complex phase cycling would be needed for more general 
states. 



III. EXPERIMENTAL METHODS 

Experiments were performed in a 7 T widebore mag- 
net with a 300 MHz Brukcr Avance Spectrometer and a 
probe tuned to 282.4MHz for 19 F measurement. A pure, 



single crystal of Fluorapatitc [Ca5(PC>4)3F] grown by flux 
method was used for the measurements [43, 44]. High 
purity of the crystal is confirmed by long relaxation time 
(Ti = 1100 s) of 19 F spins. FAp crystals have a hexagonal 
geometry with space group P63/m. The 19 F nuclei form 
linear chains along the c axis, each surrounded by 6 other 
chains. The intra-nuclear spacing within a single chain is 
d=0.3442nm and chains are separated by D=0.9367nm. 
When placed in a strong magnetic field, the nuclear spins 
interact via the secular dipolar Hamiltonian, 



3<l 



-(<jj a 1 + a^a 1 ) 

2 V x u x ' y yl 



(5) 



where the couplings depend on the relative positions as 
bji = (/io/167r)(7 2 ?i/r| ; )(l — 3cos 2 Oji), with p being the 
standard magnetic constant, 7 the gyromagnetic ratio, 
Tji the distance between nucleus j and I, and Oji the 
angle between fji and the z-axis, respectively. Due to 
1 / r 3 dependence of dipolar coupling, the largest ratio of 
in-chain to cross chain coupling is close to 40. For our 
experiments, the crystal was aligned parallel to the ex- 
ternal magnetic field in order to maximize this ratio. It 
has been shown that under these settings, and for short 
evolution times, couplings across different chains can be 
neglected and the system can be considered as a collec- 
tion of ID chains [28]. 

We performed two sets of experiments for each of the 
different initial states and readouts. First we probed the 
transport dynamics driven by the DQ Hamiltonian. For 
this purpose, the collective or end-chain magnetization 
was measured as we increased the evolution time t un- 
der the DQ Hamiltonian. We used a standard 8-pulse 
sequence [45] to implement the DQ Hamiltonian. The 
length of the 7r/2 pulse was 1.45^s. The evolution time 
was incremented by varying inter-pulse delay from 1 ps to 
6.2ps and the number of loops was increased from 1 to 
12. A recycle delay of 3000 s was used for these measure- 
ments. We restricted the evolution to a timescale where 
the ideal model applies and errors arising from leakage 
to other chains and next-nearest neighbor couplings are 
small [28]. In this timescale, the initial perturbation 
travels across « 17 spins [30], a distance much shorter 
than the average length of the chains. Thus, within the 
timescale used in this work, only the polarization starting 
from one end of the chain and moving away towards the 
other end was observed: polarization reaching the other 
end could not be observed. 

In the second set of experiments, we let the initial state 
evolve under the DQ Hamiltonian and measured the mul- 
tiple quantum coherences (MQC) to gather more infor- 
mation on the evolved state (see Appendix B for details 
on the experimental method) . The inter-pulse delay was 
varied from l^s to 6ps and the number of loops was in- 
creased from 1 to 3. We encoded coherences up to order 
4 with a K — A step phase cycling. Since these measure- 
ments involve phase cycling, a shorter recycle delay of 
1000s could be used. 
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FIG. 1: (Color online) Evolution under the DQ Hamiltonian: (a) Initial state: 8p t h. Readout: collective magnetization, (b) 
Initial state: Sp cn d- Readout: collective magnetization, (c) Initial state: Sp t h- Readout: end readout, (d) Initial state: <5p en d- 
Readout: end readout. Data points are the experimental data and lines are the fits using the analytical model described in 
appendix A. Error bars are given by the offset of the signal from zero. The fitting gives the following values for the dipolar 
coupling: 8.165 (a), 8.172(b), 8.048 (c) and 8.63 (d) xl0 3 rad/s. 



A. Experimental Results: Spin Transport 

Fig. 1(a) shows the observed evolution of the collec- 
tive magnetization under the DQ Hamiltonian starting 
from the thermal initial state (Eq. 2). The data points 
were fitted with the analytical function (Eq. Al) de- 
scribed in appendix A. In figure 1(b) we plot the sys- 
tem dynamics when starting from an end-polarized state, 
Eq. 3 (where polarization is localized at the ends of the 
chain) and reading out the collective magnetization. Fig- 
ure 1 (c) shows a complementary measurement where we 
start from thermal initial state, given by the collective 
magnetization, and read out the ends of the chains after 
evolution under the DQ Hamiltonian. Both these data 
sets were fitted by the analytical expression (A2). As it 
is evident from near perfect fitting, the analytical model 
explains the experimental data quite precisely. 

Figs. 1(a) and (b) show very different chain dynamics 
for the two initial states (with and without end selec- 
tion), giving an experimental validation of our initial- 
ization method. Furthermore, the data and fittings for 
end selection, Fig. 1(b), and end readout measurements, 
Fig. 1(c), are very similar. This indicates the robustness 
of the readout step. Finally, Fig. 1(d) shows the evo- 
lution of the end polarized initial state under the DQ- 
Hamiltonian, measured using the "end readout" scheme. 
The fitting function used is given in Eq. A3, which has 
the same form as the transport of a single excitation in a 
pure state chain [19, 34]: this experiment is thus a direct 
simulation of spin transport. 

In all the above fittings, we used the following fitting 



parameters: a scalar multiplier, a baseline constant, the 
nearest-neighbor dipolar coupling and a shift of the time 
axis. The shift of the time axis is needed since there is a 
delay of a few /is between the end of the multiple pulse se- 
quence (ideally t=0) and the measurement of the signal. 
From the fit we obtain a value for the dipolar coupling of 
8.161, 8.172, 8.048 and 8.63 x 10 3 rad/s for the four exper- 
iments, respectively. Except for the last experiment, the 
other values agree well with those obtained from similar 
measurements done on a different FAp crystal [28] and 
also with the theoretical value b=8.17 x 10 3 rad/s result- 
ing from the known structure of FAp. 

The small discrepancy in the fitting parameter for 
the last experiment, where we are initializing and read- 
ing out the chain ends, can be explained by imper- 
fections of the end-select and readout schemes. Un- 
fortunately, the phase cycling scheme does not cancel 
out zero-quantum terms. Thus, residual polarization on 
spins 2 and N-l as well as correlated states of the form 
<Jj(<Jj'_ 1 (Tj +1 + <T 7-i a j+i) contribute to errors, lowering 
the fidelity with the desired state. This effect is more 
important for the last experiment, since errors in the 
two selection steps accumulate. Moreover, the end read- 
out scheme works well only if applied to the ideal state 
which is expected after transport. The deviation of the 
prepared initial state from the ideal state further con- 
tributes to the error in the final measured data. Still, 
the agreement of the experimental data with the analyt- 
ical model indicates that these errors are small and do 
not invalidate the scheme 
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FIG. 2: (Color online) Evolution of multiple quantum coherences (0Q blue squares, 2Q red triangles, 4Q green circles): (a) 
Initial state: Spth- Readout: collective magnetization, (b) Initial state: Sp en d- Readout: collective magnetization, (c) Initial 
state: Spth Readout: end readout, (d) Initial state: 5p en d Readout: end readout. The experimental data points are fitted by 
analytical functions (blue, red and green lines) obtained from DQ Hamiltonian with NN couplings (equations B1,B2,B2 and B3 
for figures (a), (b), (c) and (d) respectively). The error bars are estimated from the deviation of 1st order quantum coherence 
from zero. Fitting of the data gives dipolar coupling: 7.971 (a), 8.077(b), 8.031 (c) and 8.492 (d) xl0 3 rad/s. 



B. Experimental Results: Multiple Quantum 
Coherences 



We present the results of the second set of experiments 
in Fig. 2, which shows the evolution of the 0, 2nd and 
4th order coherence intensities experimentally measured 
for different initial states and readouts. Fig. 2(a) shows 
the usual MQC signal, obtained measuring the collec- 
tive magnetization and starting from an initial thermal 
state (Eq. 2). The data points are fitted by analyti- 
cal functions (Bl), given in appendix B. The only vari- 
able used in these fittings was the dipolar coupling con- 
stant, for which we obtained the value b = 7.971 x 10 3 
rad/s. The data shown in Fig. 2(b) has been measured 
by first selecting the ends of the chain and then perform- 
ing the MQC measurement where collective magnetiza- 
tion is read out. Fig. 2(b), on the other hand, shows the 
data for MQC measurements starting from thermal ini- 
tial state but reading out only the spins at the chain ends. 
Both these data were fitted by the analytical functions 
(B2), giving b = 8.077 and 8.031 x 10 3 rad/s respectively. 
This is in very good agreement with the values obtained 
from the quantum transport measurements. We re- 
mark again that the results obtained from end selection, 
Fig. 2(a), and end readout measurements, Fig. 2(b), 
are very similar, thus validating the effectiveness of the 
readout step. The data for the case where we initialize 
the ends of the chains before letting the system evolve 
under DQ Hamiltonian and then read out the ends is 
shown in Fig. 2(d). Good fitting of the data with Eq. 



(B3) was obtained for b=8.492xl0 3 rad/s and shifting 
the time axis by 15/is. As mentioned earlier, we expect 
the errors in selecting the ends of the chains to add up in 
the experimental data, resulting in a slightly higher value 
of the fitting parameter b. The first two data points in 
the above mentioned figures were measured using a 4 
pulse sequence to implement DQ Hamiltonian (instead 
of a standard 8 pulse sequence), leading to greater error 
bars for these data points. 
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FIG. 3: (Color online) Logical initial state and readout of 
collective magnetization, a) Evolution under DQ Hamilto- 
nian. Data points are the experimental data and lines are 
the fits using the analytical model described in appendix 
Adequation A4. Fitting of the data gives dipolar coupling: 
7.551 xl0 3 rad/s. 
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FIG. 4: (Color online) (a) Simulated 19 F lineshapes for a chain of 11 spins, for thermal initial state (circles), ideal end-polarized 
state (triangles) and for an initial state prepared via a simulated PI sequence, with varying ti times, (b) Experimentally 
measured 19 F NMR linehsape after state initialization performed with the PI sequence for various ti time. Solid lines are for 
ti = 30^s and 35ps, which give the narrowest linewidth and the best state preparation. 



C. Quantum Information Transport 



D. 



:> F NMR lineshape analysis 



To demonstrate our ability to experimentally simulate 
not only the transport of classical information, as en- 
coded in the spin polarization, but also of quantum in- 
formation, we prepared and studied the evolution of one 
of the logical states in Eq. (4). The logical state Sp^ was 
prepared by the scheme described in section II B (simi- 
lar schemes could be used to prepare the other states in 
the operator basis). Figure 3 shows the evolution of this 
state under the DQ Hamiltonian and readout of collec- 
tive magnetization. The data points were fitted by the 
expression in Eq. (A4) giving the value of dipolar cou- 
pling constant as 7.551 xl0 3 rad/s. We note that the 
scheme for preparing this logical state involves selecting 
the ends of the chain and then creating MQC and filter- 
ing out all the terms except the double quantum terms. 
The errors involved in both these steps add up and result 
in deviation of measured data points from the analytical 
function. The experimental results, however, follow the 
expected analytical models and the dipolar coupling con- 
stants obtained from these experiments agree very well 
with those obtained from other measurements. 

In all the above described measurements where ends 
of the chains were selected and initialized, we used t\ = 
30.3^s in the pulse sequence PI. As described above, 
this value was obtained by selecting the time when po- 
larization of the spins at the ends of the chains is non 
zero while it has decayed to zero for the other spins, as 
a result of evolution under the internal dipolar Hamilto- 
nian. In order to confirm this value experimentally, we 
repeated the transport and MQC measurements at dif- 
ferent values of ti (not shown here). We observed that 
ti = 30.3^s gave the best fittings with the analytical 
functions, pointing towards the fact that the fidelity of 
end selection is highest for ti = 30.3/is. This is further 
confirmed by a detailed 19 F NMR lineshape analysis as 
described in the next section. 



A system comprising linear chains of spins, such as 
Fluor apatite, is of immense interest for NMR lineshape 
calculations by virtue of the simplicity it offers as com- 
pared to a 3-dimensional system [46, 47]. The 19 F spins in 
FAp have a characteristic 3 peak lineshape, which shows 
a strong angular dependence [48]. In our study, we uti- 
lize this angular dependence to align the crystal parallel 
to the external magnetic field in order to minimize inter- 
chain couplings. 

The measured lineshapes for thermal initial (ti = Ops) 
and end polarized (ti = 30.3/zs) using sequence P2 are 
shown in Fig. 6. Due to long Ti of 19 F nuclei, it was 



(a) 



flAA 



(b) 



-50 -30 -10 10 30 50 



30 5CH 




Frequency (kHz) 



-50 -30 -10 10 

Frequency (kHz) 




FIG. 5: (Color online) (a) The 19 F NMR lineshape (for 
ti=10/is) can be fitted by 3 Gaussian lines. Individual Gaus- 
sian lines are shown as thin dotted lines, sum of the lines is 
thick red. (b) The 19 F NMR lineshape (for ti=60/xs) can be 
fitted by 3 Gaussian lines out of which 2 lines have negative 
amplitude, (c) Amplitude of 1st Gaussian line as a function 
of ti used to fit the experimental data (red squares) and sim- 
ulated spectra (black line). 
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not possible to measure as many points as in the stan- 
dard FID measurements. This resulted in low resolution 
of the Fourier transformed signal and hence the three 
characteristic peaks of FAp lineshape are not resolved 
properly. Nonetheless, the differences in the lineshape 
and line width for the two signals are distinctly visible. 

The lineshape provides not only signatures of the sys- 
tem dimensionality, but of its initial state as well, thus we 
expect to sec qualitative differences in the NMR spectra 
for the thermal and end polarized states. We calculated 
the Free Induction Decay (FID) for a chain of N = 11 19 F 
nuclear spins evolving under the secular dipolar Hamil- 
tonian (Eq. 5) with coupling b = 8.1 x 10 3 rad/s. These 
calculations were performed also for end-polarized states, 
obtained by simulating the sequence PI for different val- 
ues of end selection time (ti). The NMR spectra ob- 
tained from these simulations are shown in Fig. 4(a), to- 
gether with the lineshapes for thermal state (Eq. 2) and 
ideal end polarized state (Eq. 3) . A few observations are 
worth pointing out. The NMR lincwidth shows a pro- 
gressive narrowing as ti is increased starting from ti = 
0/is (corresponding to the thermal state). The lincwidth 
is narrowest for ti = 30.6/is and then starts to increase 
gradually but now with an anti-phase (dispersive) com- 
ponent. A simple explanation of these features can be 
obtained by considering the signal as arising from a com- 
petition between the signal of the end-chain spins and 
the signal from the spins in the bulk. The spins at the 
chain extremities are dipolarly coupled to only one spin, 
hence we expect a splitting of the NMR line into a dou- 
blet. The spins in the chain, instead, are each coupled 
to two nearest neighbors and hence produce a broader 
NMR line, split into a triplet. As the first contribution 
increases with increasing ti time, the linewidth narrows. 
At even longer times, (ti beyond 30.6^s) two-spin cor- 
relations are created, that give rise to a dispersive spec- 
trum under the subsequent dipolar Hamiltonian evolu- 
tion. These anti-phase terms keep increasing for longer 
ti times. We observe that none of the simulated lines ex- 
actly replicates the expected lineshape for the ideal end 
polarized state. This is most likely due to less than 100% 
fidelity of end selection and initialization step. However, 
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FIG. 6: (Color online) 19 F NMR lineshape measured for ther- 
mal (black, dotted) and end polarized (red) initial state using 
pulse sequence P2. 



for ti = 30.6/iS, the width of the spectrum is very close 
to the ideal lineshape. 

The simulated lineshapes are however in good agree- 
ment with the corresponding experimentally measured 
lineshapes (Fig. 4(b)), apart from a slight asymmetry, 
which might have been introduced by a misalignment of 
the crystal with respect to the magnetic field. In particu- 
lar, the experimental spectra for ti = 30(is and 35/is have 
a width very close to the ideal end selected spectrum. 

In order to draw a quantitative comparison between 
simulated and measured NMR spectra, we fit the line- 
shapes at different ti times with a model comprising 
3 Gaussian lines, at frequencies shifted by the nearest- 
neighbor dipolar coupling. As shown in Fig. (5), both 
simulated as well as experimental lineshapes could be 
fitted reasonably well with this model. Since the outer 
lines in the Gaussian model arise only from the spins 
inside the chains, their amplitude is expected to show 
a zero crossing at ti where the chain extremities have 
maximum contribution to the measured lineshape and 
the bulk spins polarization has a very small contribu- 
tion. In our data, this is seen at 35/is indicating that end 
selection has maximum fidelity for this value of ti. This 
is in close accordance with 30.3/xs obtained from an op- 
timization of the DQ-Hamiltonian evolution fitting and 
used for our measurements. 

In the lineshape analysis described above, all the NMR 
measurements were carried out by using a 7r/2 pulse and 
measuring the resultant FID. This scheme reads out the 
collective magnetization from all the nuclei within the 
chains. It would be interesting to isolate the signal con- 
tribution from the nuclei located at the chain ends. This 
was achieved by means of pulse sequence P2, 

PI — tt/2 — (r) — tt/2 — PI, (P2) 

where the end selection step (sequence PI) is used to po- 
larize as well as read out the nuclei at the chain ends. The 
selection and readout steps are separated by a variable 
delay (r) which enables measurement of free induction 
decay (evolution under the natural dipolar Hamiltonian) 
as a function of time. The Fourier transform of this FID 
results in a NMR lineshape where only the nuclei at the 
chain ends contribute. 



IV. CONCLUSION 

In conclusion, we have studied spin transport in lin- 
ear chain of nuclear spins and experimentally demon- 
strated addressability of spins at the ends of the chain 
by means of NMR control schemes. We have shown that 
even though NMR implementation allows only collective 
control and obscrvables, we could achieve initialization as 
well as readout capabilities through a combination of co- 
herent and incoherent control. These techniques can be 
used to prepare state of relevance for quantum informa- 
tion transport as well as to monitor the dynamics of the 
end-chain spins as it evolves under the DQ Hamiltonian 
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obtained via collective manipulation of the natural dipo- 
lar interaction. We validated our method by comparing 
the experimental results with analytical solutions based 
on an idealized model, which applies to the timescales 
explored in the experiment. The good agreement of the 
data with the analytical formula confirms the preparation 
and readout of the desired state. 

In addition, we probed the states and their evolution 
by means of multiple quantum coherence measurements, 
which reveal information about multi-spin correlations. 
Again, a good conformity of the experimental results with 
the theoretical model was observed. We further opti- 
mized the end-selection scheme by a detailed analysis of 
the F 19 NMR lineshapes obtained from collective thermal 
magnetization and the end-polarized state respectively. 

Although we cannot achieve universal control of the 
end-chain spins, the initialization and readout capabili- 
ties demonstrated in this work will allow us to experimen- 
tally characterize quantum transport in spin chains. It 
will enable us, for example, to explore non-idealities that 
emerge e.g. at longer times from non-nearest neighbor 
couplings as well as couplings to other chain, and from 
the interaction of the chains with the environment. Ad- 
ditionally these methods will allow further experimental 
studies of control schemes that can enable perfect fidelity 
transfer. 
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Appendix A: Analytical solution for the evolution 
under the DQ-Hamiltonian 

Information transport in linear spin chains has been 
often studies in the literature as arising from the evo- 
lution under the isotropic XY Hamiltonian, %xy = 
J2j<e~2-(. cr j cr £ + a j a e)- For mixed-state chains, we 
showed that the DQ-Hamiltonian (Eq. 1) can as well 
drive the transport, since it is connected to XY Hamilto- 
nian by a similarity transformation. We can thus simu- 
late quantum transport with an Hamiltonian that (unlike 
the XY Hamiltonian) can be implemented experimentally 
using dipolar Hamiltonian and collective RF pulses using 
a standard sequence [36, 45]. 

The evolution of ID spin system under a DQ Hamilto- 
nian is exactly solvable in the nearest neighbor (NN) limit 
(only nearest neighbor couplings arc present and all are 
equal to b), by invoking a Jordan- Wigner mapping onto 
a system of free fermions [27, 36, 42, 45]. Various for- 
mulas describing different initial states and observables 
have been reported [20, 30, 34] and we reproduce here the 
formulas we used to interpret our experimental results in 



the main text. 

The analytical solutions for the evolution of the ther- 
mal state and end polarized state, when measuring the 
collective magnetization, are given by: 

N 

S th (t) = 5> p ,p(2*)- (Al) 
P =i 

N 

with 

oo 

= E i 2m "li 5 J2 m »+s(2bt) - i°J 2mv+a (2bt)] 

OO 

+ i 2mV [r S J2mu-5{2bt) - r a J 2m „-a{2bt)} 

Transport from one end of the chain to the other is de- 
scribed by 

S sre (t) = Aj^t) + Al N (t), (A3) 

which corresponds to the experimental situation where 
we prepared the end-polarized state and measured the 
chain ends only. Finally, we can derive the expected sig- 
nal arising from the collective magnetization when the 
initial state is the logical state 5p y L'- 

S yL {t) = A h2 {2t)+A N _ ltN (2t). (A4) 

Appendix B: Multiple Quantum NMR Spectroscopy 

Multiple Quantum spectroscopy is a powerful tech- 
nique in Nuclear Magnetic Resonance. It has the ability 
to simplify complex spectra by revealing some of the for- 
bidden transitions. Additionally, creation and evolution 
of quantum coherences can be used to probe the dynam- 
ics of a correlated many-spin system giving insight into 
dimensionality of spin system, distribution of coupling 
constants and effects of motions and quantum trans- 
port [36, 49-54]. Multiple Quantum Coherences (MQC) 
can be excited by driving the spin system by irradiation 
with cycles of multiple pulse sequences consisting of RF 
pulses and delays. For example, the so-called Double 
Quantum Hamiltonian (Eq. 1) creates even quantum co- 
herences from the longitudinal magnetization. Standard 
NMR techniques, however, measure only single quantum 
coherences. Thus, in order to probe the spin dynamics, 
it is necessary to indirectly encode the signature of the 
dynamics into single quantum coherences which can be 
directly measured inductively. This is done by labeling 
each coherence order with a different phase: If pi is the 
initial density matrix, the final density matrix p/ is given 
by: 

PS = U MQ U 4> U MQPiUl tQ UlU M Q, 
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where Umq = cxp(—iH]jQt), and = exp(— i<jxj z /2) 
is a rotation about the z axis by an angle (f>. Under this 
rotation, a coherence term of order n will acquire a phase 
nip. In order to extract the information about the dis- 
tribution of MQC, each measurement must be repeated 
while incrementing </> from to 2n in steps of S(f> = 2tt/2K 
where K is the highest order of MQC we wish to encode. 
A 7r/2 pulse is used to read out the signal at the end of 
the experiment. Finally, Fourier transform of the signal 
with respect to <fi yields the coherence order intensity: 

K 

J n (t) = Y,S k At)e- lkn5 t 

k=\ 

where S\ it) = Tr jp^)^ j i s the signal acquired in the 

fcth measurement. The technique of MQC spectroscopy 
outlined above is particularly well suited to study infor- 
mation transport by means of DQ Hamiltonian, since 
the coherence intensities show characteristics signatures 
of the occurred transport [20]. 

In a ID system with NN coupling, only zero and double 
quantum coherence orders are created by the DQ Hamil- 
tonian. Starting from the thermal initial state, the nor- 
malized intensities of the zero and DQ coherences pre- 
dicted by the analytical model are given by: 

J o /i w = ^E cos2 [ 4Mcos ^)]' 

1 k (Bl) 
4 h (t) = —^m 2 [4btco S (^ k )], 

k 

where as before N is the number of spins in the chain 
and i' k = kn/(N + 1). 

For the end-select state, zero and double quantum co- 
herence intensities given by the analytical model are as 



follows: 

J ° cnd(<) "jvTiE^ 2 ^)^ 2 ^ 008 ^)]' 

1 k (B2) 
= 1 r^- r E sin2 (^) sill2 [ 4& * C0S (^)]- 



In both Eqs. (Bl) and (B2), the normalization is chosen 
such that Jo + 2J2 = 1. 

In the case where the chain end selection is operated 
for both initialization and readout, the normalized zero 
and double quantum intensities are given by [20] : 

J ™ = JnTW ^ sin2(fc) sin2(/l) cos2[i,k + i,h) 

x (l + cos[Nfc + fc]cos[7Wi + h.]), (B3) 

Jr[t) = (jV + 1)2 E sin2 ( fc ) sin2 ( /l ) sin2 (^ + ^) 
x (1 + cos[Nk + k] cos[Nh + h]), 

where the superscript l sre' refers to select and read ends. 

Finally, for the logical state 5p y L, the MQC intensities 
are as follows: 



Jo L (t) = E sm (^) sin ( 2 V'fc) sm ( 8&tcos ^) > 

fe=l 

-1 N 

J% (t) = ——^sm(ij k )sm(2ip k ) sin (8bt cos ip k ) . 

k=l 

(B4) 
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